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We study the small-world network model, which mimics the transition between regular-lattice and 
random-lattice behavior in social networks of increasing size. We contend that the model displays a 
normal continuous phase transition with a divergent correlation length as the degree of randomness 
tends to zero. We propose a real-space renormalization group transformation for the model and 
demonstrate that the transformation is exact in the limit of large system size. We use this result to 
calculate the exact value of the single critical exponent for the system, and to derive the scaling form 
for the average number of "degrees of separation" between two nodes on the network as a function 
of the three independent variables. We confirm our results by extensive numerical simulation. 



Folk wisdom holds that there are "six degrees of sepa- 
ration" between any two human beings on the planet — 
i.e., a path of no more than six acquaintances linking 
any person to any other. While the exact number six 
may not be a very reliable estimate, it does appear that 
for most social networks quite a short chain is needed to 
connect even the most distant of the network's members 
JjJ , an observation which has important consequences for 
issues such as the spread of disease H, oscillator syn- 
chrony [||, and genetic regulatory networks Q. At first 
sight this does not seem too surprising a result; random 
networks have average vertex-vertex distances which in- 
crease as the logarithm of the number of vertices and 
which can therefore be small even in very large networks 
||. However, real social networks are far from random, 
possessing well-defined locales in which the probability 
of connection is high and very low probability of con- 
nection between two vertices chosen at random. Watts 
and Strogatz have recently proposed a model of the 
"small world" which reconciles these observations. Their 
model does indeed possess well-defined locales, with ver- 
tices falling on a regular lattice, but in addition there is a 
fixed density of random "shortcuts" on the lattice which 
can link distant vertices. Their principal finding is that 
only a very small density of such shortcuts is necessary 
to produce vertex-vertex distances comparable to those 
found on a random lattice. 

In this paper we study the model of Watts and Stro- 
gatz using the techniques of statistical physics, and show 
that it possesses a continuous phase transition in the limit 
where the density of shortcuts tends to zero. We investi- 
gate this transition using a renormalization group (RG) 
method and calculate the scaling forms and the single 
critical exponent describing the behavior of the model in 
the critical region. 

Previous studies have concentrated on the one- 
dimensional version of the small-world model, and we 
will start with this version too, although we will later 
generalize our results to higher dimensions. In one di- 
mension the model is denned on a lattice with L sites 
and periodic boundary conditions (the lattice is a ring). 
Initially each site is connected to all of its neighbors up 
to some fixed range k to make a network with average 
coordination number z = 2k M. Randomness is then in- 



troduced by independently rewiring each of the kL con- 
nections with probability p. "Rewiring" in this context 
means moving one end of the connection to a new, ran- 
domly chosen site. The behavior of the network thus 
depends on three independent parameters: L, k and p. 
In this paper we will study a slight variation on the model 
in which shortcuts are added between randomly chosen 
pairs of sites, but no connections are removed from the 
regular lattice. For sufficiently small p and large L this 
makes no difference to the mean separation between ver- 
tices of the network for k > 2. For fe = 1 it does make a 
difference, since the original small-world model is poorly 
defined in this case — there is a finite probability of a part 
of the lattice becoming disconnected from the rest and 
therefore making an infinite contribution to the average 
distance between vertices, and this makes the distance 
averaged over all networks for a given value of p also infi- 
nite. Our variation does not suffer from this problem and 
this makes the analysis significantly simpler. In Fig. |l| we 
show some examples of small- world networks. 

We consider the behavior of the model for low density p 
of shortcuts. The fundamental observable quantity that 
we measure is the shortest distance between a pair of ver- 
tices on the network, averaged both over all pairs on the 
network and over all possible realizations of the random- 
ness. This quantity, which we denote £, has two regimes 
of behavior. For systems small enough that there is much 
less than one shortcut on the lattice on average, I is dom- 
inated by the connections of the regular lattice and can 
be expected to increase linearly with system size L. As 
the lattice becomes larger with p held fixed, the aver- 
age number of shortcuts will eventually become greater 
than one and I will start to scale as log L. The transition 
between these two regimes takes place at some interme- 
diate system size L = £, and from the arguments above 
we would expect £ to take a value such that the number 
of shortcuts pk£ ~ I. In other words we expect £ to di- 
verge in the limit of small p as £ ~ p^ 1 . The quantity £ 
plays a role similar to the correlation length in an inter- 
acting system in conventional statistical physics, and its 
divergence leaves the small- world model with no charac- 
teristic length scale other than the fundamental lattice 
spacing. Thus the model possesses a continuous phase 
transition at p = 0, and, as we will see, this gives rise 
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FIG. 1. The RG transformations used in the calculations 
described in the text: (a) the transformation used for the 
k — 1 system; (b) the transformation used for k > 1, illus- 
trated in this case for k = 3. The coloring of the sites indicates 
how they are grouped under the transformations. 



to specific finite-size scaling behavior in the region close 
to the transition. Note that the transition is a one-sided 
one, since p can never take a value less than zero. In 
this respect the transition is similar to transitions seen 
in other one-dimensional systems such as ID bond or site 
percolation, or the ID Ising model. 

Barthelemy and Amaral ||] have suggested that the 
arguments above, although correct in outline, are not 
correct in detail. They contend that the length-scale £ 
diverges as 



(1) 



with r different from the value of 1 given by the scal- 
ing argument. On the basis of numerical results, they 
conjecture that r = |. Barret ||, on the other hand, 
has given a simple physical argument which directly con- 
tradicts this, indicating that r should be greater than or 
equal to 1. Amongst other things, we demonstrate in this 
paper that in fact r is exactly 1 for all values of k. 

Let us first consider the small-world model for the sim- 
plest case k = 1. As discussed above, the average dis- 
tance £ scales linearly with L for L <C £ and logarithmi- 
cally for L > If £ is the only non-trivial length-scale 
in the problem and is much larger than one (i.e., we are 
close to the phase transition), this implies that £ should 
obey a finite-size scaling law of the form 



£ = Lf(L/0, 



(2) 



where fix) is a universal scaling function with the limit- 
ing forms 



fix) ~ { 



constant 
(logx)/x 



for x <C 1 
for x > 1. 



(3) 



In fact, it is easy to show that the limiting value of fix) 
as x — > is j. A scaling law similar to this has been 
proposed previously by Barthelemy and Amaral || for 
the small- world model, although curiously they suggested 
that scaling of this type was evidence for the absence of 
a phase transition in the model, whereas we regard it as 
the appropriate form for £ in the presence of one |l(J . 

We now assume that, in the critical region, £ takes 
the form ([!]), and that we do not know the value of the 
exponent r. Then we can rewrite Eq. (^) in the form 



I = Lf(p T L), 



(4) 



where we have absorbed a multiplicative constant into 
the argument of fix), but otherwise it is the same scaling 
function as before, with the same limits, Eq. 

Now consider the real-space RG transformation on the 
k = 1 small-world model in which we block sites in adja- 
cent pairs to create a one-dimensional lattice of a half as 
many sites. (We assume that the lattice size L is even. In 
fact the transformation works fine if we block in groups of 
any size which divides L.) Two vertices are connected on 
the renormalized lattice if cither of the original vertices 
in one was connected to either of the original vertices 
in the other. This includes shortcut connections. The 
transformation is illustrated in Fig. [l]a for a lattice of 
size L = 24. 

The number of shortcuts on the lattice is conserved 
under the transformation, so the fundamental parameters 
L and p renormalize according to 
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P = 2p. 



(5) 



The transformation generates all possible configurations 
of shortcuts on the renormalized lattice with the correct 
probability, as we can easily see since the probability of 
finding a shortcut between any two sites i and j is uni- 
form, independent of i and j both before and after renor- 
malization. The geometry of the shortest path between 
any two points is unchanged under our transformation. 
However, the length of the path is, on average, halved 
along those portions of the path which run around the 
perimeter of the ring, and remains the same along the 
shortcuts. For large L and small p, the portion of the 
length along the shortcuts tends to zero and so can be 
neglected. Thus 

e' = \£ (6) 

in this limit. 

Eqs. (||) and (||) constitute the RG equations for this 
system and are exact for w> 1 and p <C 1. Substituting 
into Eq. ([|) we then find that 



r _ logjL/L') _ L 
logip'/p) 



(7) 
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Now we turn to the case of k > 1. To treat this case we 
define a slightly different RG transformation: we group 
adjacent sites in groups of k, with connections assigned 
using the same rule as before. The transformation is 
illustrated in Fig. [j]b for a lattice of size L — 24 with 
k = 3. Again the number of shortcuts in the network 
is preserved under the transformation, which gives the 
following renormalization equations for the parameters: 

l! = L/k, p = k 2 p, k' = 1, d = i. (8) 

Note that, in the limit of large L and small p, the mean 
distance I is not affected at all; the number of vertices 
along the path joining two distant sites is reduced by a 
factor k, but the number of vertices that can be traversed 
in one step is reduced by the same factor, and the two 
cancel out. For the same reasons as before, this transfor- 
mation is exact in the limit of large L and small p. 

We can use this second transformation to turn any 
network with k > 1 into a corresponding network with 
k = 1, which we can then treat using the arguments 
given before. Thus, we conclude, the exponent r = 1 for 
all values of k and, substituting from Eq. (||) into Eq. (|4|), 
the general small-world network must satisfy the scaling 
form 
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t = jf( P kL). 



(9) 



This form should be correct for L' ^> 1 and p' < 1, 
which implies that L/k 3> 1 and k 2 p <C 1. The first of 
these conditions is trivial — it merely precludes inaccura- 
cies of ±fc in the estimate of t because positions on the 
lattice are rounded off to the nearest multiple of k by 
the RG transformation. The second condition is inter- 
esting however; it is necessary to ensure that the average 
distance traveled along shortcuts in the network is small 
compared to the distance traveled around the perimeter 
of the ring. This condition tells us when we are moving 
out of the scaling regime close to the transition, which 
is governed by (0), into the regime of the true random 
network, for which (fe) is badly violated and I is known to 
scale as log L / log k [ 5| . It implies that we need to work 
with values of p which decrease as k~ 2 with increasing 
k if we wish to see clean scaling behavior, or conversely, 
that true random-network behavior should be visible in 
networks with values of p ~ k~ 2 or greater. 

We have tested our predictions by extensive numerical 
simulation of the small- world model. We have calculated 
exhaustively the minimum distance between all pairs of 
points on a variety of networks and averaged the results 
to find I. We have done this for k — 1 (coordination 
number z = 2) for systems of size L equal to a power 
of two from 128 up to 8192 and p = 1 X 10 -4 up to 
3 x 1CT 2 , and for k = 5 (z = 10) with L = 512 ... 32 768 
and p — 3 x 10~ 6 . . . 1 x 10~ 3 . Each calculation was 
averaged over 1000 realizations of the randomness. In 
Fig. U we show our results plotted as the values of Ik / L 
against pkL. Eq. (ft) predicts that when plotted in this 
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FIG. 2. Collapse of numerical data for I according to 
Eq. (^) as described in the text. Error bars are in all cases 
smaller than the data points. Note that the horizontal axis 
is logarithmic. Inset: the collapse is noticeably poorer for 



way, the results should collapse onto a single curve and, 
as the figure shows, they do indeed do this to a reasonable 
approximation. 

As mentioned above, Barthelemy and Amaral || 
also performed numerical simulations of the small-world 
model and extracted a value of r = | for the critical 
exponent. In the inset of Fig. || we show our simula- 
tion results for k = 1 plotted according to Eq. @ using 
this value for r. As the figure shows, the data collapse is 
significantly poorer in this case than for r = 1. It is inter- 
esting to ask then how Barthelemy and Amaral arrived 
at their result. It seems likely that the problem arises 
from looking at systems that are too small to show the 
true scaling behavior. In our calculations, we find good 
scaling for L/k > 60. Barthelemy and Amaral examined 
networks with k = 5, 10 and 15 (z = 10, 20, 30) so we 
should expect to find good scaling behavior for values of 
L larger than about 600. However, the systems studied 
by Barthelemy and Amaral ranged in size from about 
L = 50 to about 500 in most cases, and in no case ex- 
ceeded L = 1000. Their calculations therefore had either 
no overlap with the scaling regime, or only a small over- 
lap, and so we would not expect to find behavior typical 
of the true value of t in their results. 

It is possible to generalize the calculations presented 
here to small-world networks built on lattices of dimen- 
sion d greater than one. For simplicity we consider first 
the case k = 1. If we construct a square or (hyper) cubic 
lattice in d dimensions with linear dimension L, connec- 
tions between nearest neighbor vertices, and shortcuts 
added with a rewiring probability of p, then as before 
the average vertex-vertex distance scales linearly with L 
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FIG. 3. Collapse of numerical data for networks based on 
the square lattice in two dimensions, as described in the text. 
Error bars are in all cases smaller than the data points. 



Alternatively, we could redefine our scaling function f(x) 
so that Ik/L is given as a function oipkL d . Writing it 
in this form makes it clear that the number of vertices in 
the network at the transition from large- to small-world 
behavior diverges as (pfc)" 1 in any number of dimensions. 

Another possible generalization to k > 1 is to add con- 
nections between all sites within square or cubic regions 
of side 2k. This gives a different dependence on k in the 
scaling relation, but r still equal to l/d. 

To conclude, we have studied the small-world network 
model of Watts and Strogatz using an asymptotically ex- 
act real-space renormalization group method. We find 
that in all dimensions d the model undergoes a continu- 
ous phase transition as the density p of shortcuts tends 
to zero and that the characteristic length £ diverges ac- 
cording to £ ~ p~ T with r = l/d for all values of the 
connection range k. We have also deduced the general 
finite-size scaling law which describes the variation of the 
mean vertex-vertex separation as a function of p, k and 
the system size L. We have performed extensive numer- 
ical calculations which confirm our analytic results. 



for small L, logarithmically for large L, and the length- 
scale £ of the transition diverges according to Eq. (|l|) for 
small p. Thus the scaling form (|]) applies for general d 
also. The appropriate generalization of our RG transfor- 
mation involves grouping sites in square or cubic blocks 
of side 2, and the quantities L, p and £ then renormalizc 
according to 



Thus 



L 1 = IL, p' = 2 d p, 



log(Z/L') 



e = u. 



i 

d' 



(10) 



(11) 



logO'/p) 

As an example, we show in Fig. || numerical results for 
the d = 2 case, for L equal to a power of two from 64 up 
to 512 (i.e., a little over a quarter of a million vertices for 
the largest networks simulated) and six different values 
of p for each system size from p — 3 x 10~ 6 up to 1 x 10~ 3 . 
The results are plotted according to Eq. (Q) with r = h 
and, as the figure shows, they again collapse nicely onto 
a single curve. 

A number of generalizations are possible for k > 1. 
Perhaps the simplest is to add connections along the 
principal axes of the lattice between all vertices whose 
separation is k or less. This produces a graph with aver- 
age coordination number z = 2dk. By blocking vertices 
in square or cubic blocks of edge k, we can then trans- 
form this system into one with k = 1. The appropriate 
generalization of the RG equations (||) is then 

L' = L/k, p = k d+1 p, k' = 1, (! = I, (12) 



which gives r = l/d for all k and a scaling form of 
£=jf((pk) 1 / d L). 
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